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Abstract 

We have searched for two-body charmless hadronic decays of B mesons. 
Final states include tttt, Ktt, and KK with both charged and neutral kaons 
and pions; 717?, Kp, and K*ir; and K(f>, K*(f>, and 4><f>. The data used in this 
analysis consist of 2.6 million BB pairs produced at the T(4S) taken with the 
CLEO-II detector at the Cornell Electron Storage Ring (CESR). We measure 
the branching fraction of the sum of B° — > 7r + 7r~ and B° — ► K + ir~ to be 
(l-SloJloJ ^ 0-2) x 10~ 5 . In addition, we place upper limits on individual 
branching fractions in the rang e from 10" 4 to 10~ 6 . 

PACS numbers:13.25.Hw,14.40.Nd 



Typeset using REVT^X 



*Permanent address: University of Hawaii at Manoa 



3 



I. INTRODUCTION 



The decays of B mesons to two charmless hadrons can be described by a b — > u tree- 
level spectator diagram (Figure |l]a), or a b — > one-loop "penguin-diagram" (Figure |l]b) 
and to a lesser extent, by the color-suppressed tree (Figure |]c) or CKM-suppressed b — > dp 
penguin diagrams. Although such decays can also include contributions from b u W- 
exchange (Figure |l]d), annihilation (Figure |l|e), or vertical W loop (Figure |l]f) processes, 
these contributions are expected to be negligible in most cases. 

Decays such as B° — > 7r + 7r~ and S° — > 7r ± p T are expected to be dominated by the b — > w 
spectator transition, and measurements of their branching fractions could be used to extract 
a value for |V^,|. The decay mode B° — ► 7r + 7r~ can be used to measure CP violation in the 
B sector at both asymmetric B factories M and hadron colliders ||. Since the 7r + 7r~ final 
state is a CP eigenstate, CP violation can arise from interference between the amplitude for 
direct decay and the amplitude for the process in which the B° first mixes into a B° and then 
decays. Measurement of the time evolution of the rate asymmetry leads to a measurement of 
sin 2a, where a is one of the angles in the unitarity triangle || . If the B° 7r + 7r~ decay has a 
non-negligible contribution from the b — > dg penguin, interference between the spectator and 
penguin contributions will contaminate the measurement of CP violation via mixing [4J], an 
effect known as "penguin pollution." If this is the case, the penguin and spectator effects can 
be disentangled by also measuring the isospin-related decays B° — > 7r°7r° and B^ — > tt^it 
[J5J. Alternatively, SU(3) symmetry can be used to relate B° — > 7r + 7r~ and B° — > K + tt~ 
Penguin and spectator effects may then be disentangled || once the ratio of the two 
branching fractions and sin 2(3 |§ are measured. 

Decays such as B° — ► K + tt~ and B° K* + tc~ are expected to be dominated by the 
b —>■ sg penguin process, with a small contribution from a Cabibbo-suppressed b —>■ u 
spectator process. Interference between the penguin and spectator amplitudes can give rise 
to direct CP violation, which will manifest itself as a rate asymmetry for decays of B° 
and B° mesons, but the presence of hadronic phases complicates the extraction of the CP 
violation parameters. 

There has been discussion in recent literature about extracting the unitarity angles using 
precise time-integrated measurements of B decay rates. Gronau, Rosner, and London have 
proposed M using isospin relations and flavor SU(3) symmetry to extract, for example, the 
unitarity angle 7 by measuring the rates of B + decays to K°tt + , K + tt°, and 7r + 7r° and their 
charge conjugates. More recent publications [PHT2I1 have questioned whether electroweak 
penguin contributions (b — > 57, b — > sZ) are large enough to invalidate isospin relationships 
and whether SU(3) symmetry-breaking effects can be taken into account. If it is possible to 
extract unitarity angles from rate measurements alone, the measurements could be made at 
either symmetric or asymmetric B factories (CESR, KEK, SLAC), but will require excellent 
particle identification to distinguish between the Kn and tot modes. 

Decays such as B — > K<p and B + — > K°tt + cannot occur via a spectator process and are 
expected to be dominated by the penguin process. Measurement of these decays will give 
direct information on the strength of the penguin amplitude. 

Various extensions or alternatives to the Standard Model have been suggested. Such 
models characteristically involve hypothetical high mass particles, such as fourth generation 
quarks, leptoquarks, squarks, gluinos, charged Higgs, charginos, right-handed W's, and so 
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on. They have negligible effect on tree diagram dominated B decays, such as those involving 
b — > cW~ and b — > uW~, but can contribute significantly to loop processes like b — > sg and 
b — > dg. 

Since non-standard models can have enhanced CP violating effects relative to predictions 
based on the standard Kobayashi-Maskawa mechanism [13.13], such effects might turn out 
to be the key to the solution of the baryogenesis problem, that is, the obvious asymmetry 
in the abundance of baryons over antibaryons in the universe. Many theorists believe that 
the KM mechanism for CP violation is not sufficient to generate the observed asymmetry or 
even to maintain an initial asymmetry through cool-down JT5|] . Loop processes in B decay 
may be our most sensitive probe of physics beyond the Standard Model. 

This paper reports results on the decays B — > tttt, B — > Kit, B — > KK, B — > up, 
B — > Kp, B — >■ K*ir, B — > K(f>, B — > K*(p, and B — > 00 [JIB] . Recent observations of 
the sum of the two-body charmless hadronic decays B° — > 7r + 7r~and i^ + 7r _ [|17] and of the 



electromagnetic penguin decay -B — > i^*7 |P-8||, indicate that we have reached the sensitivity 



required to observe such decays. The size of the data set and efficiency of the CLEO detector 
allow us to place upper limits on the branching fractions in the range 10 -4 to 10~ 6 . 



II. DATA SAMPLE AND EVENT SELECTION 

The data set used in this analysis was collected with the CLEO-II detector fl9| at the 
Cornell Electron Storage Ring (CESR). It consists of 2.42 fb _1 taken at the T(4S) (on- 
resonance) and 1.17 fb _1 taken at a center of mass energy about 35 MeV below BB threshold. 
The on-resonance sample contains 2.6 million BB pairs. The below-threshold sample is used 
for continuum background estimates. 

The momenta of charged particles are measured in a tracking system consisting of a 
6-layer straw tube chamber, a 10-layer precision drift chamber, and a 51-layer main drift 
chamber, all operating inside a 1.5 T superconducting solenoid. The main drift chamber also 
provides a measurement of the specific ionization loss, dE/dx, used for particle identification. 
Photons are detected using 7800 Csl crystals, which are also inside the magnet. Muons are 
identified using proportional counters placed at various depths in the steel return yoke of the 
magnet. The excellent efficiency and resolution of the CLEO-II detector for both charged 
particles and photons are crucial in extracting signals and suppressing both continuum and 
combinatoric backgrounds. 

Charged tracks are required to pass track quality cuts based on the average hit residual 
and the impact parameters in both the r — and r — z planes. We require that charged track 
momenta be greater than 175 MeV/c to reduce low momentum combinatoric background. 

Pairs of tracks with vertices displaced from the primary interaction point are taken as K s 
candidates. The secondary vertex is required to be displaced from the primary interaction 
point by at least 1 mm for candidates with momenta less than 1 GeV/c and at least 3 mm 
for candidates with momenta greater than 1 GeV/c. We make a momentum-dependent cut 
on the 7r + 7r~ invariant mass. 

Isolated showers with energies greater than 30 MeV in the central region of the Csl 
detector, |cos#| < 0.71, where 9 is the angle with respect to the beam axis, and greater 
than 50 MeV elsewhere, are defined to be photons. Pairs of photons with an invariant mass 
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within two standard deviations of the nominal 7r° mass [^Oj are kinematically fitted with 
the mass constrained to the 7r° mass. To reduce combinatoric backgrounds we require that 
the 7T° momentum be greater than 175 MeV/c, that the lateral shapes of the showers be 
consistent with those from photons, and that | cos#*| < 0.97, where 9* is the angle between 
the direction of flight of the 7r° and the photons in the 7r° rest frame. 

We form p candidates from 7r + 7r~ or 7r + 7r° pairs with an invariant mass within 150 MeV 
of the nominal p masses. K* candidates are selected from K + tt~, K + -k°, Kg7t + or Kgir 
pairs I^TJ with an invariant mass within 75 MeV of the nominal K* masses. We form <p 
candidates from K + K~ pairs with invariant mass within 6.5 MeV of the nominal <fi mass. 

Charged particles are identified as kaons or pions according to dE/dx. We first reject 
electrons based on dE/dx and the ratio of the track momentum to the associated shower 
energy in the Csl calorimeter. We reject muons by requiring that the tracks not penetrate 
the steel absorber to a depth of five nuclear interaction lengths. We define S for a particular 
hadron hypothesis as 

dE \ _ dE I 

C dx I measured ^x I hypothesis /_. \ 

'-'hypothesis ~ \*-) 

where a is the expected resolution, which depends primarily on the number of hits used in 
the dE/dx measurement. We measure the S distribution in data for kaons and pions using 
D° — > K~7i + decays where the D° flavor is tagged using D* + — > D°tt + decays. In particular, 
we are interested in separating pions and kaons with momenta near 2.6 GeV/c. The S w 
distribution for the pion hypothesis is shown in Figure for pions and kaons with momenta 
between 2.3 and 3.0 GeV/c. At these momenta, pions and kaons are separated by 1.8 ± 0.1 
in SL. 



III. CANDIDATE SELECTION 
A. Energy Constraint 

Since the 5's are produced via e + e~ — > T(45) — > BB, where the T(4S) is at rest in 
the lab frame, the energy of either of the two 5's is given by the beam energy, E^. We 
define AE = E\ + E 2 — E^ where E\ and E 2 are the energies of the daughters of the B 
meson candidate. The AE distribution for signal peaks at AE = 0, while the background 
distribution falls linearly in AE over the region of interest. The resolution of AE is mode 
dependent and in some cases helicity angle dependent (see section III.C) because of the 
difference in energy resolution between neutral and charged pions. For modes including 
high momentum neutral pions in the final state, the AE resolution tends to be asymmetric 
because of energy loss out of the back of the Csl crystals. The AE resolutions for the modes 
in this paper, obtained from Monte Carlo simulation, are listed in Tables | and |I[ 

We check that the Monte Carlo accurately reproduces the data in two ways. First, 
the r.m.s. AE resolution for B° —>■ h + h~ (where indicates a 7r =l= or is given by 
<TAE h+h _ = v^2cr p where a p is the r.m.s. momentum resolution at p = 2.6 GeV/c. We measure 
the momentum resolution at p — 5.3 GeV/c using muon pairs and in the range p = 1.5-2.5 
GeV/c using the modes B -> ipK } B -> Dn, and B -> D*n. We find a AE + _ = 24.7±2.3±J;| 
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MeV, where the first error is statistical and the second is systematic. This result is in good 
agreement with the Monte Carlo prediction. We also test our Monte Carlo simulation in the 
modes B + -> D°tt+ and B° -> D~>k+ (where D° -> A+vr", D° -> Agvr , and D" -> i^TT") 
using an analysis similar to our B — > K*tt analysis. Again, AE resolutions for data and 
Monte Carlo are in good agreement. 

The energy constraint also helps to distinguish between modes of the same topology. 
When a real A is reconstructed as a tv, AE will peak below zero by an amount dependent 
on the particle's momentum. For example, AE for B — > K + n~, calculated assuming B 
7r + 7r~, has a distribution which is centered at —42 MeV, giving a separation of 1.7a between 
B — > K + n~ and B — > 7r + 7r~. 

B. Beam-Constrained Mass 

Since the energy of a B meson is equal to the beam energy, we use E^ instead of 
the reconstructed energy of the B candidate to calculate the beam-constrained B mass: 
Mb = J E^ — p B . The beam constraint improves the mass resolution by about an order of 
magnitude, since |p#| is only 0.3 GeV/c and the beam energy is known to much higher pre- 
cision than the measured energy of the B decay products. Mass resolutions range from 2.5 
to 3.0 MeV, where the larger resolution corresponds to decay modes with high momentum 
7r°'s. Again, we verify the accuracy of our Monte Carlo by studying fully reconstructed B 
decays. 

The M B distribution for continuum background is described by the empirical shape 

/(M B ) oc M B VT^exp (-£ (1 - x 2 )) (2) 

where x is defined as Mb/E^ and £ is a parameter to be fit. As an example, Figure § shows 
the fit for B — ► h + 7T° background from data taken below BB threshold. 

C. Helicity Angle 

The decays B — > up, B — > Kp, B — ► K*tt, and B —>■ K<p are of the form pseudoscalar 
— > vector + pseudoscalar. Therefore we expect the helicity angle, 9h, between a resonance 
daughter direction and the B direction in the resonance rest frame to have a cos 2 9h distri- 
bution. For these decays we require | cos#f/| > 0.5. 

D. D Veto 

We suppress events from the decay B + — > D°tt + (where D° — > K + tt~ or Z)° — > A^vr ) 
or 5° — > D~7i + (where D _ — > A^tt - ) by rejecting any candidate that can be interpreted as 
B — > Dtt, with a A7r invariant mass within 2a of the nominal D mass. We expect less than 
half an event background per mode from B — > Dn events after this veto. The vetoed Dtt 
signal is used as a cross-check of signal distributions and efficiencies. 
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IV. BACKGROUND SUPPRESSION USING EVENT SHAPE 



The dominant background in all modes is from continuum production, e + e~ — > qq (q = 
u, d, s, c). After the D veto, background from b — > c decays is negligible in all modes because 
final state particles in such decays have maximum momenta lower than what is required 
for the decays of interest here. We have also studied backgrounds from the rare processes 
b — > 57 and b — > utv and find these to be negligible as well. 

Since the B mesons are approximately at rest in the lab, the angles of the decay products 
of the two B decays are uncorrelated and the event looks spherical. On the other hand, 
hadrons from continuum qq production tend to display a two-jet structure. This event shape 
distinction is exploited in two ways. 

First, we calculate the angle, 9t, between the thrust axis of the B candidate and the 
thrust axis of all the remaining charged and neutral particles in the event. The distribution 
of cos 9t is strongly peaked near ±1 for qq events and is nearly flat for BB events. Figure |] 
compares the cos 9t distributions for Monte Carlo signal events and background data. We 
require |cos#t| < 0.7 which removes more than 90% of the continuum background with 
approximately 65% efficiency for signal events ||22j| . 

Second, we characterize the event shape by dividing the space around the candidate 
thrust axis into nine polar angle intervals of 10° each, illustrated in Figure the i th interval 
covers angles with respect to the candidate thrust axis from (i — 1) x 10° to i x 10°. We fold 
the event such that the forward and backward intervals are combined. We then define the 
momentum flow, X{ (i = 1,9), into the i th interval as the scalar sum of the momenta of all 
charged tracks and neutral showers pointing in that interval. The 10° binning was chosen to 
enhance the distinction between BB and continuum background events. 

Angular momentum conservation considerations provide additional distinction between 
BB and continuum qq events. In qq events, the direction of the candidate thrust axis, 9 q q, 
with respect to the beam axis in the lab frame tends to maintain the 1 + cos 2 9 qq distribution 
of the primary quarks. The direction of the candidate thrust axis for BB events is random. 
The candidate B direction, 9b, with respect to the beam axis exhibits a sin 2 9b distribution 
for BB events and is random for qq events. 



A Fisher discriminant |23| is formed from these eleven variables: the nine momentum 
flow variables, | cos 6^1, and | cos^l- The discriminant, J 7 , is the linear combination 
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T = s £ j a i x i (3) 

i=i 

of the input variables, Xj, that maximizes the separation between signal and background. 
The Fisher discriminant parameters, «j, are given by 

<*i = f:K + U* j )- 1 x(tf-^. (4) 
J=l 

where Ufj and are the covariance matrices of the input variables for signal and background 
events, and fij, ^ are the mean values of the input variables. We calculate using Monte 
Carlo samples of signal and background events in the mode B — > 7r + 7r~. 
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Figure |6| shows the T distributions for Monte Carlo signal in the mode B° — > 7r + 7r~, 
and data signal in the modes B — > Dtt. Figure ^ also shows the T distributions for Monte 
Carlo background in the mode B — > h + TT~ and below-threshold background data for modes 
comprising three charged tracks or two charged tracks and a 7r°. The T distribution for signal 
is fit by a Gaussian distribution, while the T distribution for background data is best fit by 
the sum of two Gaussians with the same mean but different variances and normalizations. 
The separation between signal and background means is approximately 1.3 times the signal 
width. We find that the Fisher coefficients calculated for B° — > 7r + 7r~ work equally well for 
all other decay modes presented in this paper. Figure [7] shows the remarkable consistency 
of the means and widths of the T distributions for signal and background Monte Carlo for 
the modes in this study. 

V. ANALYSIS 

For the decay modes B — > txtx , B — > Kit, and B — > KK, we extract the signal yield using 
a maximum likelihood fit. For the other decay modes, we use a simple counting analysis. 
Both techniques are described below. 

A. Maximum Likelihood Fit 

We perform unbinned maximum likelihood fits using AE, Mb, T , and dE/dx (where 
appropriate) as input information for each candidate event to determine the signal yields for 
B° -> 7T+7T-, K+iT-, K+K~, 7r°7r°, K n , and B + -> tt+tt , K+n°, K°ir + . Five different 
fits are performed as listed in Table ||. 

For each fit a likelihood function £ is defined as: 

N 

C = n P (A, f m ; (A£, M B , T , dE/dx),) (5) 

8=1 

where P (fi, f m ; (AE, Mb, T, dE/dx),) is the probability density function evaluated at the 
measured point (AE, Mb, T , dE/dx)i for a single candidate event, i, for some assumption 
of the values of the yield fractions, fj, that are determined by the fit. iV is the total number 
of events that are fit. The fit includes all the candidate events that pass the selection criteria 
discussed above as well as | cos 6^1 < 0.7, and < T < 1. The AE and M B fit ranges are 
given in Table 0. 

For the case of B — > h + h~, the probability Pi = P (f±, f m ; (AE, Mb, T, dE/dx)/) is 
then defined by: 

Pi = J"' 7T7T P^TT + fKirPK-K + IrrPrK + ~ fnTT ~ fl<TT ~ f K k) P (^) 

pC _ f C pC , f C pC , /i _ fC _ fC )p c 
r — J-K-K r -K-K + J Kir Kir + I 1 J tttt J K-w) r KK 

where, for example, P^ n (P^ n ) is the product of the individual probability density functions 
for AE, Mb, J 7 , and dE/dx for 7r + 7r _ signal (continuum background). The signal yield in 
B° — * 7T + 7r _ , for example, is then given by N n7T = f^ n x N. 
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The central values of the individual signal yields from the fits are given in Table [HI 
None of the individual modes shows a statistically compelling signal. To illustrate the fits, 
Figure |8| shows Mb projections for events in a signal region defined by \AE\ < 2a a_e and 
T < 0.5 and Figure |9| shows the AE projections for events within a 2a Mb cut and T < 0.5. 
The modes are sorted by dE/dx according to the most likely hypothesis and are shown in the 
plots with different shadings. Overlaid on these plots are the projections of the fit function 
integrated over the remaining variables within these cuts. (Note that these curves are not 
fits to these particular histograms.) 

Our previous publication |U| reported a significant signal in the sum of B° — > 7r + 7r _ and 
B° — > K + 7i~. While our current analysis confirms this result, we now focus on separating 
the two modes. We separate the systematic errors that affect the total yield from those 
that affect the separation of the two modes. We do this by repeating the likelihood fit using 
N sum = N n7T + N Kn , R = N wn /N sum , and fixing N KK = 0, its most likely value. We find: 

AT — 17 9+5.6 +2.2 
Jv sum — J- 1 4.9 —2.5 

R = 0.54±g;|g ± 0.05 

where the first error is statistical and the second is systematic (described below). The result 
of this fit is shown in Figure [H| This figure shows a contour plot (statistical errors only) of 
N sum vs. R in which the solid curves represent the na contours (n =1-4) corresponding to 
decreases in the log likelihood by 0.5n 2 . The dashed curve represents the 1.28a contour, from 
which estimates of the 90% confidence level limits can be obtained. The central value of N SVLm 
has a statistical significance of 5.2a. The significance is reduced to 4.2a if all parameters 
defining £ are varied coherently so as to minimize N sum . Further support for the statistical 
significance of the result is obtained by using Monte Carlo to draw 10000 sample experiments, 
each with the same number of events as in the data fit region but no signal events. We then 
fit each of these sample experiments to determine N sum in the same way as done for data. 
We find that none of the 10000 sample experiments leads to N snm > 10. 

None of the physical range of R can be excluded at the 3a level. However the systematic 
error of R is only 10% (see below and Table [TV] ) . We therefore conclude that our analysis 
technique has sufficient power to distinguish the 7t + tt~ mode from K + tt~, but at this time 
we do not have the statistics to do so. 

Since none of our fits has a statistically significant signal, we calculate the 90% confidence 
level upper limit yield from the fit, N 90 , given by 

I N9 °C^(N)dN 

f Q °° £ max (N)dN U ' yU [i) 

where C max (N) is the maximal C at fixed N to conservatively account for possible correlations 
among the free parameters in the fit. The upper limit yield is then increased by the systematic 
error determined by varying the parameters defining C within their systematic uncertainty 



as discussed below. Table [TTJ summarizes upper limits on the yields for the various decay 
modes. 

To determine the systematic effects on the yield due to uncertainty of the shapes used 
in the likelihood fits, we vary the parameters that define the likelihood functions. The 



variations of the yields are given in Table IV. The largest contribution to the systematic 
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error arises from the variation of the Mb background shape. For this shape, /(Mg) oc 
— x 2 exp(— £(1 — x 2 )) (x = M B /E h ), we vary E h by ±1 MeV, consistent with observed 
variation; we vary £ by the amount allowed by a fit to background data (below-threshold 
and on-resonance AE sideband) which pass all other selection criteria. To be conservative, 
we allow for correlated variations of E\> and £. 



B. Event-Counting Analyses 

In the event- counting analyses we make cuts on AE, Mb, T , and dE/dx. The cuts for 
AE and Mb are mode dependent and are listed in Table |. We require JF < 0.5. Tracks 
are identified as kaons and/or pions if their specific ionization loss, dE/dx, is within three 
standard deviations of the expected value. For certain topologies, candidates can have 
multiple interpretations under different particle hypotheses. In these cases we use a strict 
identification scheme where a track is positively identified as a kaon or a pion depending on 
which dE/dx hypothesis is more likely: we sort the modes with two charged tracks plus a 7r° 
(7r + p~, 7r°p°, K + p~ , K* + tt~, and i^*°7r°) by requiring strict identification for both charged 
tracks. For modes with three charged tracks (vr + p°, K + p°, and K* 0, k + ) we require strict 
identification of the two like-sign tracks, while the unlike-sign track P3| is required to be 
consistent with the pion hypothesis within two standard deviations. We separate modes 
with one charged track plus two 7r°'s (p + 7r° and K* + tt°) by requiring strict identification of 
the charged track. 



Figures |TT]-|14] show M B distributions for B — > up, B — > Kp, B — > K*ir, B — > K(f>, B — > 
K*(p and B — > (jxp candidates (after making the cuts on AE, T , and particle identification 
described above.) The numbers of events in the signal regions are listed in Table [V|. 

In order to estimate the background in our signal box, we look in a large sideband region 
in the AE vs. M B plane: 5.20 < M B < 5.27 GeV and \AE\ < 200 MeV. The expected 
background in the signal region is obtained by scaling the number of events seen in the 
on-resonance and below-threshold sideband regions (weighted appropriately for luminosity). 
Scale factors are found using a continuum Monte Carlo sample which is about five times the 
continuum data on-resonance. In many modes, the backgrounds are so low that there are 
insufficient statistics in the Monte Carlo to adequately determine a scale factor. For these 
modes, we calculate upper limits assuming all observed events are signal candidates. The 
estimated background for each mode is also listed in Table 0. 

Although we find that there are slight excesses above expected background in some modes, 
no excess is statistically compelling. We therefore calculate upper limits on the numbers of 
signal events using the procedure outlined in the Review of Particle Properties PDA for 



evaluation of upper limits in the presence of background. To account for the uncertainties in 
the estimated continuum background we reduce the background estimate by its uncertainty 
prior to calculating the upper limit on the signal yield. 



VI. EFFICIENCIES 

The reconstruction efficiencies were determined using events generated with a GEANT- 
based Monte Carlo simulation program p5l . Systematic uncertainties were determined using 
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data wherever possible. Some of the largest systematic errors come from uncertainties in the 
efficiency of the | cos9 T \ < 0.7 cut (6%), the uncertainty in the 7r° efficiency (7% per 7r°), and 
the uncertainty in the Kg efficiency (8% per Kg). In higher multiplicity modes, substantial 
contributions come from the uncertainty in the tracking efficiency (2% per track). In the 
B — > 7rp, Kp, K*tt analyses, the simulation of the efficiency for the particle identification 
method has a systematic error of 15%. For the event-counting analyses, the uncertainty in 
the T < 0.5 cut is 5%. 

The total detection efficiency, £ , is given by £ = £ r x £ d , where £ r is the reconstruc- 
tion efficiency and £ d is the product of the appropriate daughter branching fractions. The 
efficiencies, with systematic errors, are listed in Tables |TTT| and [VT| . 

VII. UPPER LIMIT BRANCHING FRACTIONS 

Upper limits on the branching fractions are given by N\jl/(£Nb) where iV\jL is the upper 
limit on the signal yield, £ is the total detection efficiency, and Nb is the number of B 0, s 
or B +, s produced, 2.6 million, assuming equal production of charged and neutral B mesons. 
To conservatively account for the systematic uncertainty in our efficiency, we reduce the 
efficiency by one standard deviation. The upper limits on the branching fractions appear in 
Tables [TTI| and 0. 

VIII. SUMMARY AND CONCLUSIONS 

We have searched for rare hadronic B decays in many modes and find a signal only 
in the sum of tt + it~ and K + tt~ . The combined branching fraction, B(tt + tt~ + K + n~) = 
(l.Slg^loJiO^) x 10~ 5 , is consistent with our previously published result. We have presented 
new upper limits on the branching fractions for a variety of charmless hadronic decays of 
B mesons in the range 10~ 4 to 10~ 6 . These results are significant improvements over those 
previously published. Our sensitivity is at the level of Standard Model predictions for the 
modes ir + ir~, K + n~, tt + tt°, K + n°, 7r ± p =F , K + <fi, and K*°(J). 
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TABLES 



TABLE I. Resolutions of AE and the signal regions for AE and AM B = M B - 5280 MeV for 
the event-counting analyses. Indicated in parentheses are the K* decay modes used. 











Mode 




\AE\ 


IAMrI 




(MeV) 


(MeV) 


(MeV) 


lt^p T 


25-46 


< 2<r a 


< 6.0 


n p 


46 


<- on 

\ i7U 


\ u.u 


71 P 


AO 


<r 50 

\ <JU 


/fin 

\ U.U 


7T p 


ou 


\ 1UU 


\ D.U 


J\ p 


AO t-tu 




\ u.u 


j\ p 


99 

A A 


\ UU 


\ u.u 


n. p 


93 

AO 


\ tJU 


\ u.u 


j\ p 


A A *-to 




\ u.u 










(K+TT ) 


25-40 


< 2a a 


< 6.0 


(K°TT+) 


21 


< 50 


< 6.0 


K*°TT° 








{1\ TT ) 


44 




\ u.u 


I\ Ji 








(K + TT°) 


50 


< 100 


< 6.0 


(K°TT + ) 


45 


< 90 


< 6.0 


I\ n 








[J\ IT ) 


93 

AO 


<r 50 

\ tJU 


/fin 

\ U.U 


<K°TT°) 


22-40 




< 6 


K°(j) 


18 


< 45 


< 6.5 


K+cp 


23 


< 60 


< 6.5 


K*°(/) 








(K+TT-) 


20 


< 50 


< 6.5 




24 


< 60 


< 6.5 


K*+4> 








(i^+vr ) 


23 


< 60 


< 6.5 




17 


< 45 


< 6.5 




16 


< 40 


< 6.5 




a The AE resolution and cut 


are functions of the helicity angle. 
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TABLE II. Resolutions of AE, the fit regions in AE and Mb, and the number of events, N, in 
the fit regions for the likelihood analyses. 

Fit region 





&AE 


AE 


M B 




Mode(s) 


(MeV) 


(MeV) 


(GeV) 


N 


ir + ir-/K + iT~/K+K- 


±25 


-185 < AE < 140 


5.21 < M B < 5.30 


453 




+43/-55 


±300 


5.20 < M B < 5.30 


896 


7r°7r° 


+51/-85 


±300 


5.20 < M B < 5.30 


104 


K°TT° 


+44/-53 


±200 


5.20 < M B < 5.30 


44 


K°TT+ 


±27 


±200 


5.20 < M B < 5.30 


220 



TABLE III. Results from the likelihood analyses: the signal yield central value from the fit 
(Ns), detection efficiencies (£), the 90% confidence level upper limit on the number of signal events 
(-ZVtjl), the 90% CL upper limit of the branching fraction (UL B), and the theoretical predictions 
for the branching fractions [p^-28|. We also include the measured branching fraction (B) for the 



sum of 7r + 7r and K + tt , where the first error is statistical, the second is the systematic error from 
the yield, and the third is the systematic error from the efficiency. 







£ 


B 




UL B 


Theory 


Mode 


N s 


(%) 


(10" 5 ) 




(io- 5 ) 


(IO" 5 ) 


h+TT' 


1 7 r) + 5.fc> + 2.2 

i '- z -4.9-2.5 


37 ±3 


1.81-1^ ± 0.2 








7T + 7r~ 


9.411? 


37 ±3 




17.9 


2.0 


1.0-2.6 


K+-K- 


7.9±1;I 


37 ±3 




15.3 


1.7 


1.0-2.0 


K+K- 


o o+ - 8 
u - u -o.o 


37 ±3 




3.5 


0.40 




K+TT° 


4.91U 


33 ±3 




11.2 


1.4 


0.3-1.3 


vr+vr 




33 ±3 




13.1 


1.7 


0.6-2.1 


7T 7r° 


1 2 +L7 


26 ±4 




5.2 


0.91 


0.03-0.10 


K°TT + 


c o+3.5 
d - z -2.8 


11 ± 2 




11.3 


4.8 


1.1-1.2 


K°TT° 


O q+2.2 
z -' 5 -1.5 


7± 1 




6.2 


4.0 


0.5-0.8 
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TABLE IV. Dominant variations in the upper limit signal yield (%) due to systematic uncer- 
tainties in the fit shapes. 



Mode 


Background Mb 


Signal M B 


Signal AE 


T 


dE/dx 


Total 




8.6 


2.9 


5.4 


4.7 


5.1 


13 


K+it- 


5.0 


3.3 


2.6 


2.5 


4.6 


8 


K+K- 


4.7 


< 1 


< 1 


2.2 


3.2 


9 


J\ 7T 


o.y 


< u.o 


Z.O 


z.U 


O 1 

z. / 




7T+7T 


13.5 


2.8 


3.3 


3.2 


1.7 


15 




12.1 


3.1 


1.2 


5.1 




14 


K°ir+ 


7.5 


4.8 


1.8 


5.7 




9 




6.9 


1.0 


1.1 


1.6 




7 


a N 

1 'sum 


+10.4/-12.8 


+3.9/-3.5 


+1.5/-0.7 


+5.9/-6.4 


±1 


+13/-15 


a R 


+2.0/-3.3 


+1.7/-1.4 


±7.2 


+1.7/-2.6 


±5.6 


+9/-10 


a Systematic errors on central value. 
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TABLE V. Results of the event-counting analyses: the number of events in the signal region 
(Ns), the estimated background in the signal region (JYb), the 90% confidence level upper limit on 
the branching fractions (UL B), and theoretical predictions [26-28]. Indicated in parenthesis are 
the K* decay modes used. 



Mode 


N s 


N B 


UL B 

(io- 5 ) 


Theory 
(IO" 5 ) 


7T ± /J =F 


7 


2.9 ±0.7 


8.8 


1.9-8.8 


7T O° 


1 


1.8 ± 0.6 


2.4 


07-0 23 


7T+P 
r 


4 


2.3 ± 0.3 


4.3 


0.0-1.4 


TrV 
r 


8 


5.5 ± 1.2 


7.7 


1.5-3.9 


K + p~ 


2 


2.0 ± 0.4 


3.5 


0.0-0.2 


K°n° 


n 


n 


3 9 


004-0 04 


K+p° 


l 


3.8 ± 0.2 


1.9 


0.01-0.06 


K°p+ 

n r 


o 





4.8 


0-0.03 




3 


0.7 ± 0.2 


7.2 


0.1-1.9 




(3) 


(0.7 ±0.2) 






(K°7r+) 


(0) 


(0) 








o 


1.1 ± 0.3 


2.8 


0.3-0.5 




(0) 


(1.1 ± 0.3) 








4 


1.9 ± 0.7 


9.9 


05-0 9 




(3) 


(1.9 ± 0.7) 






(K°TT + ) 


(1) 


(0) 






K*°ir + 


2 


1.0 ± 0.6 


4.1 


0.6-0.9 




(2) 


(1.0 ± 0.6) 






(K°ir°) 


(0) 


(0) 






K°cp 


1 





8.8 


0.07-1.3 


K + cp 








1.2 


0.07-1.5 




2 





4.3 


0.02-3.1 


(K+TT-) 


(2) 


(0) 








(0) 


(0) 






K* + (f> 


1 





7.0 


0.02-3.1 


(K+TT°) 


(0) 


(0) 






(K°7T + ) 


(1) 


(0) 






66 








3.9 
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TABLE VI. Reconstruction efficiencies (£ r ), the products of the appropriate daughter branching 
fractions (£d), and total detection efficiencies {£ = £ r x £d) for the event-counting analyses. 



Mode 

7r°p° 
vr+p 

7T°p+ 



K+p- 
K°p° 
K+p° 
K°p+ 

K* + TT~ 

(K+TT°) 

(K°7T+) 
K*°7T 

(K+TT-) 
K* + 7T° 

(K+TT°) 

(K°vr+) 

(K+TT-) 
(K°7T ) 

K°4> 
K + (p 

(K+7T-) 
(K°7T ) 

K*+(f> 

(K + TT°) 



£ r (%) 

5.3 ±1.1 
6.5 ± 1.3 

7.4 ± 1.5 

5.5 ± 1.1 

5.7 ± 1.1 

7.8 ± 1.2 
7.1 ± 1.4 
6.4 ± 1.0 



4.5 ±0.8 

9.8 ±2.0 

6.1 ± 1.2 

3.9 ±0.8 

7.6 ± 1.5 

7.1 ± 1.4 
7.9 ± 1.6 

11.9 ± 1.8 
17.8 ±2.7 

16.2 ±2.4 
5.6 ±0.8 

9.2 ± 1.4 
9.4 ± 1.4 
11.0 ± 1.7 



0.988 
0.988 
1.0 
0.976 



5.2 ± 1.0 
6.4 ±1.2 
7.4 ±1.5 
5.4 ± 1.1 



0.988 
0.343 
1.0 
0.339 



5.7 ±1.1 
2.7 ±0.4 

7.1 ± 1.4 

2.2 ±0.3 



0.329 
0.228 



3.7 ±0.4 
(1.5 ±0.2) 
(2.2 ±0.3) 



0.657 

0.325 
0.226 

0.665 
0.113 



4.0 ±0.8 

3.0 ±0.4 
(1.3 ±0.2) 
(1.7 ±0.3) 

5.6 ±0.9 
(4.7 ±0.9) 
(0.9 ±0.2) 



0.168 
0.491 

0.327 
0.055 

0.162 
0.112 
0.241 



2.0 ±0.3 

8.7 ±1.3 

5.6 ±0.8 
(5.3 ±0.8) 
(0.3 ±0.1) 

2.6 ±0.5 
(1.5 ±0.2) 
(1.1 ±0.2) 

2.7 ±0.4 
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FIG. 1. Feynman diagrams for rare hadronic B decays: (a) b — > u external W emission, (b) 
b — > s,d loop or gluonic penguin, (c) b — > u internal W emission, (d) b — > u W exchange, (e) 
annihilation, and (f) vertical W loop. 



20 



120 I 1 r 



1 DD - 



0990995-017 



t 1 1 r 



6D - 



4D - 



Pions 



■ ■ Kaons 



nf 



_l I U I I L 



FIG. 2. Distribution of S n for kinematically identified high momentum kaons and pions from 
D* + — ► D°7T + ; D° — ► K~7r + decays. The solid line shows SV for pions and the dashed line shows 
S n for kaons. 
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FIG. 3. Mb distribution from below-threshold background events (squares) and the fit to the 
parameterization given in the text (curve). The mass for the below-threshold data is shifted up to 
match the kinematic endpoint of the on-resonance data. 
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FIG. 4. The cos 9t distributions for background data (squares) and B° — > tt + tt Monte Carlo 
signal (histogram). 
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FIG. 5. Illustration of the first three of the nine polar angle intervals. 
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1.DD 



FIG. 6. The T distribution for B° — > 7r + 7r~ Monte Carlo (solid histogram), B — > Z)7r signal 
data (filled squares), the fit to the signal data (solid curve), the background Monte Carlo (dotted 
histogram), background data (open squares), and the fit to the background data (dotted curve). 
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FIG. 7. The means of the Fisher output distributions for signal (filled circles) and background 
(open circles) for the modes in this study. The error bars indicate the width of the distributions. 
Since the backgrounds in the <j) modes are small, their background means are poorly measured. 
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FIG. 8. M B plots for (a) B° -► vr+vr- (unshaded) 5° -► K+vr- (grey), and 5° -► K+K~, 
(black) (b) 5+ -► vr+vr (unshaded) and 5+ -► K+tt° (grey), (c) 5° -► 7r°7r°, (d) B° -► ^°7r°, 
and e) .B — > i^°7r + . The projection of the total likelihood fit (solid curve) and the continuum 
background component (dotted curve) are overlaid. 
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FIG. 9. AE plots for (a) B° -► vr+vr- (unshaded) B° -► K+vr- (grey), (b) 5+ -► vr+vr 
(unshaded) and 5+ -► K + ir° (grey), (c) £° -► 7r°7r°, (d) B° -► K°7r°, and (e) 5 -► £T 7r + . In plots 
(a) and (b), the projection of the total likelihood fit (upper solid curve), the irir signal component 
(dotted curve), the Kir signal component (dashed curve), and the background component (lower 
solid curve) are overlaid. In plots (c)-(e), the total likelihood fit (solid curve) and background 
component (dotted curve) are overlaid. 
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FIG. 10. The central value (+) of the likelihood fit to N sum = N n7V + Nk-k and R = N^/Ng^ 
for B° — > 7r + 7r _ and B° — > K + 7r~. The solid curves are the n<r contours, and the dotted curve is 
the 1.28(7 contour. 
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FIG. 11. M B plots for (a) B° -► vr ± p^, (b) B° -► ttV, (c) -► and (d) B + -► vr p + - 

The signal region is indicated. 
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FIG. 12. M B plots for (a) B° -► K+p~, (b) B° -► K°p°, (c) 5+ -► and (d) B+ -► 

The signal region is indicated. 
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FIG. 13. M B plots for (a) B° -► Jf+Ti--, (b) 5° K*°7r°, (c) £+ -► K*+tt°, and (d) 
5 + — > ivT* 7r + . The shaded events are from X* — ► K^n decay modes and the unshaded events are 
from K* — > i^ + 7r decay modes. The signal region is indicated. 
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FIG. 14. M B plots for (a) B° -► (b) 5+ -► K + <f>, (c) 5° -► iT*V, (d) B + -► ir"+</>, and 

(e) 5° — > 00. The shaded events are from K* — > -K^tt decay modes and the unshaded events are 
from if* — > K + tt decay modes. The signal region is indicated. 
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